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Abstract 

A simple analytical solution is found to the Dirac equation for the combina- 
tion of a Coulomb potential with a linear confining potential. An appropriate 
linear combination of Lorentz scalar and vector linear potentials, with the 
scalar part dominating, can be chosen to give a simple Dirac wave function. 
The binding energy depends only on the Coulomb strength and is not affected 
by the linear potential. The method works for the ground state, or for the 
lowest state with I = j — \, for any j. 



Typeset using REVTgX 



♦Internet address: V5030E@VM.TEMPLE.EDU 



1 




It is often useful in quark model calculations to have a simple wave function to help 
in parameterizing physical quantities. A variety of non-relativistic and relativistic wave 
functions have been used for this purpose, but the desire for simplicity has often limited the 
physical reasonableness of the wave function. In this paper, we propose a simple relativistic 
wave function that is the solution of the Dirac equation for a potential that has the general 
features that have been suggested by QCD. The potential we use has a Coulombic part and 
a linear confining part. The following calculation is exceedingly simple, but we have not 
seen it presented anywhere. We present it here so that it might be useful to anyone who 
would like to use a simple, but reasonable, relativistic potential and wave function. 

The Dirac equation we solve is 

[a • p + j3m — — + P/j,r + vr]ip = Eip, (1) 
r 

where ex. and (3 are the usual Dirac matrices. The linear confining part has been split into 
two parts, Pfir which is a Lorentz scalar potential, and ur, a Lorentz vector potential. We 
will see that the inequality /i > v is necessary for the potential to lead to a bound state. 
That is, the Lorentz scalar part must dominate over the Lorentz vector part of the confining 
potential. The four component wave function ip can be written in terms of two component 
spinors u and v as 

(2) 

with iV an appropriate normalization constant. The two component spinors u and v satisfy 
the equations 

(<T-p)u = [E + m + - + (fi- v)r\v (3) 
r 

(a ■ p)v = [E — m-\ (/! + u)r]u. (4) 

For u, we choose the function 

u = r h - x e- ar e--^ r \, (5) 
where x is a constant two component spinor. This is just the form of the exact Coulomb 
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solution with an added factor e~z a r . Acting on this form of u with the operator cr -p gives 
the result that 

(a ■ p)u = %{a ■ r)[a + + a 2 r]u. (6) 

Comparision of Eqs. (3) and (6) for (cr ■ p)u shows that v can be written as 

v = ry(<7 • r)u, (7) 
where 7 is a constant factor, given by 
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a 1 — b a 2 

7 = = = • (8) 

m + E A \i- v v 7 

The three equalities in Eq. (8) result from comparing the constant term, the - term, and 
the r term in Eqs. (3) and (6) separately. The form for v given by Eq. (7) can now be 
substituted into Eq. (4). This results in the same form for u as in Eq. (5), and leads to three 
new equations for the constant 7 



m — E A /i + v 



The relations in Eqs. (8) and (9) can be combined, resulting in 

b = VT^ (10) 

a = m\ (11) 
E = my/1 - A 2 . (12) 

These are just the same results, including the same energy, as for a simple Coulomb potential 
without the added linear confinement. The confining potential constants must be related by 



v 



-fiVT^x 2 , (13) 



so that the vector part of the linear potential must be less in magnitude than the scalar 
part. The Gaussian constant a 2 is given by 

a 2 = A/i, (14) 

so that fi must positive and v negative. Note that the scalar confining strength can take 
any value, as long as ji is positive and the two constants, \i and v, are related by Eq. (13). 
We see that the solution is not completely general, but can only work if the three potential 
constants satisfy the constraint of Eq. (13). However this still permits a wide range of 
confining potentials, because the sum // + v is unconstrained. Only the ground state will 
have the simple analytic form of Eqs. (5) and (7). The radially excited states could not be 
written so simply, and their energies would depend on the confinement constants \i and v 
as well as on A . 

The method also works for the lowest excited state for which I — j — \- In this case, 
equations (11)-(14) get modified by the substitution 

A -> A/k, (15) 

where k = j + \ is the principal quantum number of the state. This can be seen by writing 
the Dirac wave function in terms of radial functions f(r) and g(r) 

/ f(r)y j 1 \ 

^ = n[ M / 5* , (16) 

where y\ is a two component angular spinor function corresponding to total angular mo- 
mentum j and orbital momentum I. The radial functions satisfy the equations 
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U— V=- 

ar r I 



E + m H h (/J. — f)r 



' d 1 + k\ 
dr r j ^ 



A 



E — m H (// + f) 7 " 

r 



For the function f, we choose the form 



/ = r 6 - 1 e _ar e-5 ar2 . 
Then substituting this form into Eqs. (17) and (18) leads to the result 

9 = if 

with 



k — 6 A 



a 



7 = 



m-E 



A + 6 ij, — v a 2 m + E 

Equation (21) leads to the result 

b = V/t 2 - A 2 , 



(17) 
(18) 

(19) 
(20) 

(21) 
(22) 



and the modification A — > A//c in Eqs. (11)-(14). 

In summary, we have presented a simple Dirac wave function for a Coulomb potential 
with linear confinement. 
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